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Abstract 

Finite time coherent sets [ ] have recently been defined by a measure based objective function describing the 
degree that sets hold together, along with a Frobenius-Perron transfer operator method to produce optimally 
coherent sets. Here we present an extension to generalize the concept to hierarchically define relatively coherent 
sets based on adjusting the finite time coherent sets to use relative measure restricted to sets which are developed 
iteratively and hierarchically in a tree of partitions. Several examples help clarify the meaning and expectation 
of the techniques, as they are the nonautonomous double gyre, the standard map, an idealized stratospheric flow, 
and empirical data from the Mexico Gulf during the 2010 oil spill. Also for sake of analysis of computational 
complexity, we include an appendix concerning the computational complexity of developing the Ulam-Galerkin 
matrix estimates of the Frobenius-Perron operator centrally used here. 

Keywords: Coherent Pairs, Relative Measures, Relatively Coherent Structures, Frobenius-Perron 
Operator, Subdivision Flow Chart. 

1 Introduction 

Central to understanding mixing and transport mechanisms is the related question of defining partitions relative 
to which the transport can be discussed. To this end, the concept of almost invariant sets defined for autonomous 
systems [4, 2, 7] and coherent sets for non autonomous systems [ ] are central since transport may be defined as 
the measure of the set that leaves the partition element corresponding to the almost invariant set (or finite time 
coherent set) in a given time epoch. See also [ ]. Transfer operator methods are proven to be computationally 
effective for use in identifying almost invariant sets and finite time coherent sets. See [ , , ]. 

The method to identify coherent pairs used here is based on the Frobenius-Perron operator, the Ulam-Galerkin 
method and the thresholding method. See [ ]. 

In this work, we extend the concept of finite time coherent pairs to incorporate relative measure, and we call this 
relatively coherent pairs. This extension provides the theoretical framework to simply apply the definition of finite 
time coherence iteratively and at each stage we hierarchically define the relative measure on each element of the 
sub partitions which are developed. The results can be collected in a tree structure to emphasize the hierarchical 
nested nature of such partitions. 

Also we present a adapted thresholding method that can identify the finite-time relatively coherent structures in 
successive scales of time-dependent systems. We study four examples in this paper, for both closed systems and open 
systems. In the first and second example, we show our method identifies the fine scaled relatively coherent structures 
in a nonautonomous double gyres and a standard map. In the third example, an idealized zonal stratospheric fiow, 
our method gives fine scaled details with respect to relative coherence. In our fourth example, we show the method 
is also efficient at an open system, corresponding to oceanographic ffows in the Mexico Gulf. 

This paper is organized as follows. In Sec. 2, we define relative coherent structures. Then we briefiy describe 
the Frobenius-Perron operator, and the Ulam-Galerkin matrix estimate and the thresholding method. In Sec. 3, we 
give the details of the algorithm and a successive fiow chart for explanation. From Sec. 4 to Sec. 7, we apply the 
method to four examples. Conclusions are given in Sec. 8. In Sec. 9, the appendix, we analyze the computational 
complexity. 

2 Relatively coherent structures 

We define relatively coherent structures by specializing the definition of coherent pairs, which we now review 
Let (1], A, /i) be a measure space, where ^ is a a-algebra and /i is a normalized measure that is not necessarily 
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invariant. The key to the speciahzation to relatively coherent pairs is the use of an iteratively defined relative 
measure, refined from the initial globe measure /i. 

Generally, we assume that Q C M^. Given a time-dependent fiow <l>(z,t;r) : 1] x R x R ^ 1], through the time 
epoch r of an initial point z at time t. A coherent pairs (A^, A^+r) can be considered as a pair of subsets of ft such 
that. 

Definition 2.1 [ is a {po^t^r) -coherent pair if 

p^{At, At^r) := K^t n ^{At^r^t + r, -r))//i(At) > po (2.1) 

where the pair (At^At-^T-) are ^robust ^ to small perturbation and p{At) = p{At-^T-). 

Note that the definition centrally depends on the full measure /i on cj, and we will substitute successive relative 
measures on refinements. Now we consider a relative measure on u induced by /i, where a; is a nonempty measurable 
subset of ft. In this way enters refinements of the initial partition on successive scales, a relative measure of w 
to Q is 

""•-'-^ (2.2, 

for all Ae A. 

From the above definition, it follows that the space (cj, A\uj^ Pu) is also a measure space, where A\uj is the restric- 
tion of cj to ^ and p^ is a normalized measure on uj. We call the space {uo^ A\oj^ Poj)^ the relative measure space. 
Now, we define the relatively coherent pairs. 

Definition 2.2 Relatively coherent structures are those {po^t^r) -coherent pairs defined Definition 2.1, with respect 
to given relative measures on a subset uo (ZVt, of a given scale. 

To find relatively coherent structures in time-dependent dynamical systems, the basic tool is the Frobenius- 
Perron operator. Let (Vt^A^p) be a measure space and /i is a normalized Lebesgue measure. If : 1] ^ 1] is 
a nonsingular transformation such that p{S~^{A)) = for all A G ^ satisfying p{A) = 0, the unique operator 
P : L\Q) L^Q) defined by, 

Pf{x)p{dx) = / f{x)p{dx) (2.3) 

for all A G ^ is called the Frobenius- Perron operator corresponding to S, where f{x) G L^{Q) is a probability 
density function. See In our case, S can be considered as the fiow map <l> and the formula above can be written 
as 

Pt,rf{z) := f{S-\z)) • \det D{S-\z))\ = /($(^, t + r; -r)) • \det D^z, t + r; -r)|. (2.4) 

Suppose X is a subset of M, let Y" be a set that includes S{X). We develop partitions for X and Y respectively. In 
other words, let {5^}^^ be a partition for X and {Cj}^^^ be a partition for Y. The Ulam-Galerkin matrix follows 
a well-known finite-rank approximation of the Frobenius-Perron operator, which is of the form 



_ n{B,ns-\c,)) 

where /i is the normahzed Lebesgue measure on $7. As usual, we numericahy approximate Aij by, 

p _ #{xfc -.XkeBjk Sjxk) e Cj} 

where the sequence {x^} is a set of test points (passive tracers). See [6]. 



The following thresholding method from [ ] finds optimally coherent pairs in a time-dependent dynamical sys- 
tem, with respect to the chosen measure. We will iteratively adapt the method to the relative measure. This 
algorithm thresholds to the singular values and singular vectors of the matrix P obtained by the Ulam-Galerkin 
method : 
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i. Calculate the second singular value and corresponding left and right singular vectors {xi} and {yj} of the 
Ulam-Galerkin matrix. 



ii. Find values {{bk^Ck)} as pairs such that, 

p{X{b,), Y{c,)) = ^"^'g-^^-^-^^-^'^^ (2.7) 

by thresholding, where 

Pi = fi{Bi),X{bk) = yJi:x.>hkBi and Y{cu) = U,:^^.>c,C,. (2.8) 

iii. Choose a partition related to a pair (6*,c*) of (bk^Ck) such that 

p* =max{ p{X{bk),Y{ck))}. (2.9) 

k 

The partition is maximally coherent with respect to /i on 1] and the test set {(6/e,c/e)}. 



3 Algorithm 

We now describe how to find relatively coherent pairs. By the thresholding Eqs. 2.7 - 2.9, we have obtained optimal 
coherent pairs, which are defined as (Xi,Yi) and (X2,l2)- Yi can be considered as the image of Xi under a fiow 
^ in time-r. In order to find relatively coherent structures in Xi,yi, X2 and ¥2^ we define relative measures on 
each of these sets. Define relative measures fiXi{S) and vyi{T) on the coherent pair (Xi,yi), according to Eq. 2.2. 
Then we have the measure spaces (Xi, , /^Xi) and {Yi^ A\yi^vyi) with jj^Xi and vyi the normalized measures ^. 

The relative measures for both Xi and Yi allow the adaptation of the thresholding methods on (Xi, , MXi) 
and {Yi^ A\yi^vy^) under the relative measures. Then follows two relatively coherent pairs in the previous coherent 
pair (Xi,Fi), which are now named as (Xii,Yii) and (Xi2,yi2). Also, the coherent pair (X2,l2) can be divided 
to two relatively coherent pairs, (X2i,l2i) and (X22,^22)- 

Next, we repeat the building process above, that is, to define eight normalized measures on each of Xn, Yn, X12, 
Y12, X21, I21, ^22 and I22, respectively, such that they become new spaces with corresponding relative measures, 
we then apply the adapted thresholding method on these new spaces to get more relatively coherent structures. 

Now, we state our hierarchical adaptation method as an algorithm, and for convenience, we use (X^,y^) and 
(Xj^Yj) to denote two coherent pairs, i and j which can be stated to emphasize the hierarchy tree. 

Algorithm 1 

1. Define relative measures juxi and 'Uy. and relative measure spaces (X^, A\xi^ MxJ and (F^, ^|y., vyJ for (X^, Y^), 
where 



f^xAS) = TTTTT and vy,{T) = ——- (3.1) 

for all S cXi and for all T C Yi. 

2. Apply the adapted thresholding method on {Xi^ A\xi^ fJ^xJ and (F^, ^|y. , 'Uy.) for (X^,y^) to develop refined 
coherent pairs {Xa^Ya) and (Xij^Yij). 

3. Repeat the above two steps for the pair (Xj^Yj) to obtain further refined coherent pairs (Xji^Yji) and 

^The probability measure v can be considered as the discretized image of 12. The detail of construction of v can be found in [8] 
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According to the first three steps, a relatively coherent structure can be denoted by 



(3.2) 



and kp G = 1,2, ...q, through g-steps of the algorithm, q -\- 1 levels deep into the tree. See Figure 1. 

That is, the subscript kik2...kq can be any possible finite q permutations of i and j. We usually choose i = 1 
and j = 2. Figure 1 is a fiow chart depicting four levels, which describes the steps to find relatively coherent 
pairs in finer 'scales'. The chart emphasizes a 'tree' structure. 

However, we cannot repeat the algorithm forever, so we must decide a stopping criterion. The following step 
is as a completion for the algorithm. 

4. Stop a given branch if in Eq. 2.9, 



where po G (0, 1) is a threshold from Eq. 2.1, descriptive of optimal coherence which is not very coherent. 
The Gulf Example in Figure 7 on how such stopping criterion leads to the number of coherent pairs less than 2^. 



p* =max{ p{Xk^k2...k,{h)^ykik2...k,{ck))}' 

k 



(3.3) 



is such that 



p < Po- 



rn 
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Figure 1: Successive relatively coherent sets tree as per Algorithm 1. 



4 Example 1. The Nonautonomous Double Gyre. 

Consider the nonautonomous double gyre system, 

X = — 7rAsin(7r/(x, t)) cos(7r?/) 

y = 7TAcos{7Tf{x^t)) sm{7Ty) (4-1) 

where f{x^t) = esm{ujt)x^ + (1 — 2esm{ujt))x^ e = 0.25, = 27r and A = 0.25. See [^]. 

Let the initial time be t = and the final time be t = 10. We use 24,200 identical triangles ^ to cover 

the region [0, 2] x [0, 1] by Delaunay triangulation. Then randomly and uniformly we choose 10,000,000 points in 
the region as our initial conditions. For good sampling, the relationships between the given grid and the necessary 
number of points we choose is discussed in Appendix I, on computational complexity. We numerically calculate 
the final status of these points by the Runge-Kutta method to estimate the flow. Because the double gyre model 
is an area-preserving system, the same triangulation can be used as the image partition {Cj}'j^i^. The Ulam- 
Galerkin's transition matrix estimates the Frobenius- Perron operator that has the size 24200 by 24200. We apply 
the thresholding method on the matrix to find the first two coherent pairs in the initial status and final status. See 
Figure 2. The first two coherent pairs are colored blue and red in the first level of both the upper and lower charts. 
We define the left part of the initial status as Xi and the left part of the final status as Yi, both of which are filled 
with blue as halves of the initial status and final status. Thus, two relative measures can be defined on each of 
these two parts separately. 
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(4.2) 



On the other hand, we can do the same in the right red-fihed regions which we cah the initial status X2 and the 
final status Y2 to develop another two relative measures as follow, 



^^^^^) = M^' '^'^^^^^ = W2)' ^'-'^ 

Now Xi and Yi can be considered as the initial status and the final status of a refined relative sub- 'system'. 
Likewise for X2 and Y2. By the same process as with the whole double gyres system, we can get some new coherent 
structures in the new system consisting of Xi and Yi. In Figure 2, following the first blue arrow between first level 
and second level of both fiow charts, Xi, the blue half on the first level of the upper chart is divided by blue and 
light blue, we define the blue part as Xn and the light blue part as X12 in the second level. Correspondingly, we 
have Yii which is blue and Y12 which is light blue in the second level of the lower chart of Figure 2 as the outcome 
states of Xii and X12. Xn and Yn are relatively coherent structures, and so are X12 and Yi2- 

As above, we can develop an X21 that is red and an X22 that is light green from X2 in the second level of the 
upper chart of Figure 2; and Y21 that is red and Y22 that is light green from Y2 in the second level of the lower 
chart of Figure 2. The same subscript means X21 and Y21 are a relatively coherent pair, so are X22 and Y22' Now 
we have four relatively coherent pairs in the second level. 

We can eventually get the tree structures in Figure 2 by repetition of the process. There are 8 relatively 
coherent structures shown in the third level with different colors and 16 relatively coherent structures in the fourth 
level with different colors. In Figure 2, we can see the egg shaped relatively coherent structures which are four 
resonance 'islands' as expected in such hamiltanian twist maps, [ ]. Even finer structures will be revealed by further 
refinement and sufficient sampling to allow appropriate resolution. Appropriate sampling in a given refinement scale 
is discussed in terms of computational complexity in developing a given Ulam-Galer kin's matrix for a given fine 
grid. 
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2"^ level 




3^^ level 4^^^ level 

(a) t=0 






l^t level 



2"^ level 



level 



4^*^ level 



(b) t=10 



Figure 2: The upper figure is the initial status of the double gyre and the lower figure is the image under the time-r 
flow, where r = 10. The same color areas between the two associates one relatively coherent pair. In this case, we 
have 16 different relative coherent structures with different colors. By following the colored arrows, we can see the 
relative coherent structures through four levels of reflnement. 
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5 Example 2. The Standard Map. 



Consider the standard map, 

Pn+l =Pn ^Ksm{On) 

On+1 = On +Pn+1 (5.1) 

where Pn and On are taken modulo 27r. See [ ], and this is a map on the torus Q = [0, 1) x [0, 1). We wih study 
the case that K = 1.2 in our example, as this is well known to he K = 1.2 > Kcr = 0.971635 shortly after the 
breakup of the last 'golden' torus allowing momentum boosting orbits and a mixed chaotic and ordered phase space 
including periodic elliptic islands, [ ]. 

We choose 10,000,000 points, randomly and uniformly in the area [0,1) x [0,1). We set X := [0,1) x [0,1) 
and Y := [0, 1) x [0, 1) and then uniformly triangulate to cover both the initial status and the mapping image. 
Therefore, enumerating the triangulation grid {Bi}'^^-^^ and {Cj}^2=T^5 initial points do 10 iterations, after 

which we obtain the final status of the system and the Ulam-Galer kin's matrix is 20,000 by 20,000. 

The first level of Figure 3 shows two large coherent pairs in the Standard Map, which are filled by blue and red, 
respectively. By prior knowledge on this benchmark problem, it is clear that the boundary between the primary 
blue and red coherence estimates the cantorus remnant of the golden area resonance, an expected act came as this 
is known to remain a primary pseudo-barrier to transport when K = 1.2 still not much larger than Kcr = 0.971635. 

Iteratively repeating the process according to the algorithm for each of Xi, Yi, X2 and Y2 yields 4 relatively 
coherent structures which are colored blue, light blue, red and green in the second level. The third level in Figure 
3 tells us there are 8 different such structures in total in the first three levels. The outcome partition here shows 
a familiar depiction of the resonance layers known to be due to cantorus pseudo-barriers which cause the famously 
slow transport for the standard map, [15]. 
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level level 3^^ igygi 

(a) At the beginning 




V' level level 3^^ ^^^^^ 

(b) After 10 iterations 

Figure 3: Standard Map Eq. 5.1, structured relative coherence hierarchy tree arranged as in (a) - (b). 
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6 Example 3. An Idealized Stratospheric Flow. 



Next, consider the Hamiltonian system 



where 



dx/dt = -d^/dy 
dy/dt = d^/dx, 



^(x, y^ t) = C2,y — UoLtanh{y /L) + A^JJ^Lsech^ {y / L)cos{kix) 

-\-A2UoLsech'^ {y / L)cos{k2X — (J2t) 
-\-AiUoLsech'^ {y / L)cos{kix — ait) 



(6.1) 



(6.2) 



This is a quasiperiodic system that represents an ideahzed zonal stratospheric flow, [ ]. There are two known 
Rossby wave regimes in this system. Let Uq = 63.66, C2 = 0.205/7o, c3 = 0.7/7o, A3 = 0.2, A2 = 0.4, Al = 0.075 and 
the other parameters in Eq. 6.2 be the same as stated in [ ]. 

We choose 20,000,000 points in the domain X = [0, 6.37l7r * 10^] x [-2.5 * 10^, 2.5 * 10^] of the flow and 
use 32,640 triangles as the partition {Bi}^^^^ for the initial status points and 39,694 triangles as the partition 
{Cj}^^i^ for the final status of the points. Note that this system is 'open' relative to the windows X chosen. The 
two coherent pairs are colored blue and red which are defined as (Xi^Yi) and (X2,l2) in the first level of Figure 
4. Again, we now build the relative measures and tree of relatively coherent pairs. By applying the method as we 
have done with the previous two examples, we develop four and eight different coherent structures for the second 
level and the third level, respectively. In Figure 5, we repeat the color scheme in both the initial status and final 
status in each of the second and third levels for the relative coherent structure. Thus, we now see a much finer 
scaled relative coherence in the dynamical system than previously seen. 





(a) t=0 



(b) t=0 



M 10 





(c) t=10 days 



(d) t=10 days 



Figure 4: Relative coherence in the Rossby system Eq. 6.2. The fist level partition (left hand side) and third 
level partition (right hand side) of both the initial and final status of the zonal flow. Compare to the hierarchical 
structure as emphasized in Figure 5. 
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7 Example 4. The Gulf of Mexico. 



In our last example, we consider the Mexico Gulf. The data is the same as was used in [ ] and formed by the 
method in [ ] and [ ] . The difference between the Gulf model and the above three examples is the Gulf is a open 
system, that is, there is water entering and exiting the region. See Figure 6. This is the reason why in Figure 7, 
the shapes of the whole Gulf water of the initial status and final status are slightly different at the bottom region 
and top region. 




-98 -96 -94 -92 -90 -88 -86 -84 -82 -80 -78 



(a) 

Figure 6: Vector field describing surface flow in the Gulf of Mexico on May 24, 2010, computed using the HYCOM 
model [HYCOM, 2010]. Note the coherence of the Gulf Stream at this time. Oil spilling from south of Louisiana 
could flow directly into the Gulf Stream and out towards the Atlantic. This is an open system relative to this 
window shown. Horizontal and vertical units are degrees longitude (negative indicates west longitude) and degrees 
latitude (positive indicates north latitude), respectively. 

We choose 20,000,000 points uniformly and randomly in the water region as the initial status, on more details 
of data. See [ ]. The final status is the positions of these points after 6 days. We use 21,645 triangles as a 

partition of X and 21,210 triangles {Cj}^'^^^ as a partition of Y . After applying our subdivision method on these 
triangles, the results are shown in Figure 7. In this example, we set the po = 0.9998 as the threshold the stopping 
criterion. Therefore, the number of relatively coherent pairs is not equal 2^, where q is defined in Eq. 3.2. 
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p* = 0.9997 < po 



p* = 0.9999 > Po 



p* = 0.9999 > Po 
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p* = 0.9992 < Po 
^0 
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(a) t=0 

p* = 0.9997 < Po 
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:^HP p* = 0.9999 >po V ^ 
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1 



p* = 0.99978 < Po 

yjp*=0 _ . 



= 0.99989 > Po 



p* = 0.9999 > Po 



4. 



1st level level 3^^ level 4^^^ level 

(b) t=6 days 

Figure 7: Hierarchical relative coherence in the Gulf of Mexico following the flow according to vector field data as 
illustrated in Figure 6. Tree structure and relatively coherent pairs coloring as in (a)-(b). 
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8 Conclusions 



We have defined a concept of relative coiierence based on relative measure, as a generalization of coherent pairs. 
We also have introduced a recursive method of detecting relatively coherent structures under flows in a finite time, 
which is based on relative measures, with respect to the restricted Frobenius- Perron operator. Relative measures 
are used to build a hierarchy of relatively coherent pairs in successive levels, which can be illustrated in a natural 
tree structure of relative coherence. 

We have demonstrated the method with the double gyres, the standard map, a Rossby wave system and data 
from the Gulf of Mexico. These examples have included hierarchical structure, open and closed systems, a system 
known only through data, and use of the stopping criterion. 



9 Appendix I: On Computational Complexity. 

In this section, we analyze the computational complexity to properly develop an Ulam-Galerkin matrix. The 
question is how many initial points should we use for a given grid. In other words, how many initial points is 
enough to well represent the whole domain for a given grid. Intuitively, we may wish to add 'as many as possible'. 
Generally, finer grids require exponentially more points, depending on the dimensionality, but also on the local 
stretching of the map. The following discussion is base on the Lipschitz constant and Gronwall's inequality. See 
[16]. There are some other works related to this topic. See [ , 9, 10, 13]. 

Consider a grid that consists of square boxes with length q. The triangulation we use is handled similarly, but 
rectangles will simplify this discussion even if triangles allow for the powerful Delaunay triangulation algorithms in 
practice. For convenience, we use boxes instead of triangles, but one box can be easily changed to two triangles by 
cutting through the diagonal. 

There are several steps to find how many points we need. Let f : X x [to^t] ^ X be a flow and X be compact. 

I. To find an e such that for Vxi,X2 G X, \xi{to) — X2{to)\ < e we have 

\xi{t) - X2{t)\ <q, (9.1) 

where xi(to) and ^2(^0) are the initial status of xi,X2 and xi(t),X2(t) are the final status of xi,X2. The 
purpose of this step is to avoid that the distance d between the final positions of two closed initial points under 
the flow change dramatically, in our case, it means d > q. See Figure 8 case (a). So, we must control e such 
that case (b) holds for all initial points. 

By Gronwall's inequality, 

\xi{t) - X2{t)\ < \x,{to) - X2(to)|e'^l*-*°l, (9.2) 
where M is the Lipschitz constant. Assuming / is uniformly continuously different able. 



M = max 

{x,t)eXx[to,t] 



(9.3) 



We must control the distance between xi{t) and X2{t) through controlling the distance between xi{to) and 
X2{to). Let 

(9.4) 



gM|t-to| * 

Then we have, 

\xi{t) - x{t)\ < \xi{to) - x(to)|e^l^-^°l < £6^1^-^°! = q, (9.5) 



for all X G X satisfying \xi{to) — x{to)\ < £. That is, for any x{to) in an e-ball of xi{to), x{to) will keep its 
distance with Xi{to) less than q through time t — to. 

II. Next, we consider the more general case of the e-ball of xi{to); consider e-balls centered on each sample data 
X which must cover the whole X. Since X is compact, the e-balls of a finite subset {xk}^=i cover X. Let 
{xk}k=i be initial conditions of the flow, with images under the flow {xk{t)}^^i in a finite time t — to. Define, 

/ = min {distance between each of four boundaries of some j box and Xk(t)} 
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The Initial Status The Final Status 

(a) d > q. A failure scenario 



The Initial Status The Final Status 

(b) d < q. The success scenario 



Figure 8: (a)Two points in box i evolve to boxes j + 1 and j — 1 and without further samples, the resulting Ulam- 
Galerkin matrix would miss a transition i ^ j which is expected in a continuous transformation /. (b) The success 
scenario we wish to guarantee. Points xi and X2 are close enough that the / never casts them further than across 
successive rectangles j — I and j. Recall that a continuous / will map a connected set to a connected set and this 
should be reflected in the images of the coarse representation of boxes cast across boxes. 




(a) some j box from X(t) 

Figure 9: The red arrow indicates the minimum distance between a point Xk{t) for some k and four boundaries of 
a box j for some j that contains Xk{t). We must check all in {xk{t)}^^-^ with all the boxes of the partition of 
X{t) to develop 

where the j box is as defined above, is the image of X which we denote X{t). See Figure 9. 
Here, again, we use Gronwall's inequality as follows, 

\x(t) - Xk{t)\ < \x{to) - Xfe(to)|e^l*-*°l < h 

where k = 1,2, ... P. Let 

I 



oM\t-to\ ' 



(9.6) 



we have 



\xk{t) - x(t)\ < \xk{to) - a;(to)|e^l*-*°l < ee^l*"*"! = I, 



(9.7) 
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which is to say, ah points x in an ^''-ball of Xk in X for some k wih map into the /-bah of Xk{t) in X(t) for 
the same k. 

III. In this step, we build a new cover for X{t) with /-bahs by the fohowing. There are two types of bahs of the 
cover. Type I covers are those /-bahs never touch the boundaries of the grid; Type II covers are those /-bahs 
centered at boundaries and their centers never in the /-bahs of Type I, see Figure 10. 




TYPE I 



TYPE II 



(a) two types of covers 



Figure 10 



X{t) is compact, so we have a finite subset of this new cover stih covering X(t). Let the centers of these bahs 
be {xq(t)}^^^^ for Q is some integer. Since the flow is assumed continuous as is its inverse is continuous, (see 
[16]), the pre-image of the these {xq{t)^ / j-bahs is stih a cover of the pre-image X of which is also finite. 

By the above discussion, more specifically, the finite cover of X consists of ^''-balls of some of {^g(^)}^=i- That 
is, the pre-image of Type I /-balls in X{t) and another type of covers whose image are Type II /-balls in 
the shape of these covers may be vary. Thus, in X, we have the following case, see Figure 11. 

Pre-image of TYPE I 



Pre-image of TYPE II 

of covers 

Figure 11: We use an ellipse to denote the pre-image of some Type II /-ball, but the pre-image may be any shape 
which depends on the inverse of the fiow. 

Consider the areas which are not overlapping of the balls of both of X and X(t). The fiow is a bijection 
between these two areas. See Figure 12, which is to say, the points in red ^''-balls in X will never go out of the 
/-balls in X(t). 

Notice the important connection between the non-overlapping areas of X and we can control the non- 

overlapping area in X by shrink the size of /-balls in X(t). In other words, if / become smaller, so is the 
non-overlapping area in X, that is, the non-overalpping area in X can be ignored for some small /. Then 
we can guarantee a sufficient sampling for a given grid with length q if we choose mesh-grid of points with 




(a) the pre-image of the two types 
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distance e' or smaller. The number of points we need in each box of the partition of the initial status X should 
be more than ^ . To obtain a e'', we just choose a small / which depends on the Lipschitz constant of a specific 
problem. Note that this description is in terms of a uniform e' grid. However, a uniformly random cover where 
this e' condition is satisfied is also sufficient. 




(a) the relationships between the non-overlapping areas 

Figure 12 
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